Signal Analysis based on Complex Wavelet Signs 

L. Demaret 1 , P. Massopust 1 ' 2 -*, M. Storath 1 ' 2 



Abstract 

We propose a new analysis tool for signals that is based on complex wavelet signs, called a signature. The 
complex- valued signature of a signal at some spatial location is defined as the fine-scale limit of the signs of its 
complex wavelet coefficients. We show that the signature equals zero at sufficiently regular points of a signal 
whereas at salient features, such as jumps or cusps, it is non-zero. At such feature points, the orientation of 
the signature in the complex plane can be interpreted as an indicator of local symmetry and antisymmetry. 
We establish that signature is invariant under fractional differentiation and rotates in the complex plane 
under fractional Hilbert transforms. Thus, the signature may be regarded as a complementary object to the 
local Sobolev regularity index. We derive an appropriate discretization, which shows that wavelet signatures 
can be computed explicitly. This allows an immediate application to signal analysis. 

Keywords: Wavelet signature, Complex wavelets, Signal analysis, Hilbert transform, Fractional 
differentiation and integration, Phase, Sobolev regularity, Singular support, Randomized wavelet 
coefficients, Rademacher sequence, Salient feature 
2010 MSC: 42C40, 94A12, 44Af5 



1. Introduction 

The determination and discrimination of salient features, such as jumps or cusps, is an important task 
in signal processing. Classical approaches assume the interesting features of a signal to be points of low 
regularity. In this context, local regularity is measured in terms of the (fractional) differentiability order, 
e.g., in the sense of local Holder, Sobolev or Besov regularity. Since such measures of smoothness only rely 
on the modulus of wavelet coefficients, they do not take into account wavelet sign (or phase) information. 

We may observe the shortcomings of a purely modulus based approach by considering the two functions 
f(x) = sgnx and g(x) = 2 log |a:|. Since / and g are related by the Hilbert transform, their wavelet coefficients 
are equal with respect to the order of magnitude. Hence, the locally symmetric singularity of / and the 
locally antisymmetric singularity of g at the origin cannot be distinguished using a purely modulus-based 
signal analysis. 

We present a new signal analysis tool, which is based on the sign of the wavelet coefficients. To this end, 
we investigate the fine scale limits of the signs of the wavelet coefficients 

af{b) := lim sgn (/, n a>b ) := lim \{' Ka ' b \ , (1) 

where k is a complex-valued wavelet, a > the scale, and be 1 the location. The complex-valued quantity 
crf(b) is called the signature of f at location b. We shall see that the signature allows the local analysis 
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of isolated salient features. Hereby, the orientation of the signature within the complex plane may be 
interpreted as an indicator of local symmetry or antisymmetry. In particular, we show that the signature is 
purely imaginary at a jump, whereas it is purely real at a cusp. Moreover, the signature is invariant under 
fractional Laplacians, i.e., 

a((-A)tf) = af, 

and it serves as a multiplier when acting on the fractional Hilbert transform, in the sense that 

a(n a f)=e^af. 

Therefore, the signature may be interpreted as being "dual" to the local Sobolev regularity index, which is 
invariant under fractional Hilbert transforms but shifts under fractional Laplacians. We also establish that 

sing supp / <£_ supp af and supp af <£ sing supp /. 

Thus, a singularity in the classical sense need not coincide with a signature- type singularity. Finally, the 
signature classifies each point of a signal, whose signs are randomly disturbed, as a non-salient point. 

We also propose a method to numerically compute the signature of real world signals and validate the 
theoretically developed concepts by numerical experiments. 

1.1. Related approaches and our contribution 

The local decay rate of the modulus of wavelet coefficients is often used to characterize the singular 
support of a function. It even provides a means of extracting the local regularity. (See, for instance, 
[1].) The local decay rate is also used in the context of detection and classification of salient features in 
signals and images. Attempts along these lines can be found in the book by Mallat, see [2, Chapter 6.2-6.3] 
and references therein, where the procedure is called wavelet transform modulus maxima and applied to 
multiscale edge detection. However, modulus based methods do not, and can not, take into account the 
information contained in the sign (or phase) of the wavelet coefficients. 

This stands in stark contrast to the importance of sign information in the signal structure. For instance, 
it was shown in [3] that bandpass signals are characterized, up to a constant, by their zero crossings, which 
in turn are determined by the changes in the sign of the function values. Thus, when dealing with bandpass 
wavelets, the wavelet coefficients are determined by their sign. In this context, the well known experiments 
by Oppenheim and Lim [4] indicate that sign information is more important for the reconstruction of images 
than the modulus. 

The sign information of the Fourier coefficients is used in signal analysis via phase congruency [5, 6]. 
In [7], Kovesi added a multiscale aspect to phase congruency and applied it successfully to signal analysis. 
However, the approach was not derived rigorously and involves several heuristic methods and principles. 
Regarding the signs of complex wavelets, it has been observed in [8] and [9] that the lines of constant phase 
in the continuous wavelet transform of a Gaussian complex wavelet converge towards the singularities in the 
signal. However, to the authors' knowledge, these observations have not been cast into rigorous statements. 

We introduce a new approach to incorporate the sign of complex wavelet coefficients into signal analysis. 
In particular, we use the complex-valued signature to delineate singularities in signals. This novel method 
is based on complex wavelets, which are complex-valued functions whose real and imaginary parts are 
Hilbert transforms of each other. Complex wavelets have been successfully applied in various signal and 
image processing applications, for instance in denoising [10, 11]. In this paper, we use the class of complex 
bandpass wavelets that have a one-sided and positive Fourier spectrum (Definition 1). This class includes, 
for instance, the complex Meyer-type wavelets. The main idea behind the signature of a signal / is to 
investigate the fine scale behavior of the complex wavelet signs, sgn (/, K a ,b), defined as 

It \ (/) K a,b) 
Sgn /, K a , b ) = — 

I if, K a.b) | 

if (/, K a ,b) 7^ and zero otherwise. The signature af at a location b £ R is equal to the limit 

af(b) = lim sgn (/, n a . b ) , 

a— >0 



if the limit exists for all wavelets k in our class and if the limit is independent of the choice of k; otherwise, 
we set <jf(b) = 0; see Definition 2. We consider b € K to be a salient point of a signal / if the signature at 
b is not equal to 0. Indeed, if / has a jump discontinuity at b, we prove in Theorem 12 that the limit equals 
±i, whereas if / is of cusp-type at b, we observe in Example 9 that the limit equals ±1. Furthermore, if 
the signal is locally polynomial around b, i.e. regular, then either the above limit does not exist, or equals 
zero. (Theorem 5 and Corollary 6) These results show that signature is capable of distinguishing between 
different types of singularities and may be used to classify them. Hereby, the orientation of the signature 
within the complex plane can be seen as an indicator of local symmetry or antisymmetry. This coincides 
with the analogous interpretations for phase congruency [5, 6, 7]. A similar connection between the fine 
scale wavelet phase and the local geometry of a singularity was observed in [12, Chapter 4.3]. 

We set the signature into relation to the classical local smoothness analysis in terms of Sobolev regularity. 
The local Sobolev regularity index of a function is invariant under the application of the fractional Hilbert 
transform but changes under the action of the fractional Laplacian [13]. We show that the signature is "dual" 
to the local Sobolev regularity index in the sense that the fractional Hilbert transform acts like a multiplier 
on erf, whereas of remains invariant under the application of the fractional Laplacian; see Theorem 14 and 
Theorem 15. 

In [14] it is shown that random modifications of the wavelet signs significantly alter the shape of a signal. 
Since the moduli of the wavelet coefficients remain unchanged so does the classical modulus based regularity. 
In Theorem 19 we show that such a random modification of the signs of the wavelet coefficients yields a 
signature identically equal to zero. Hence, a signal of randomly disturbed wavelet signs does not possess 
any salient features with respect to our definition of signature. 

In order to compute the signature of a real world signal /, we introduce a discretization based on the 
mean values 

1 N 1 N (f « ) 

^ := ivg sgn ^^> = ivgK7^- (2) 

where {o,j}°°^ 1 is a discrete zero sequence of scale samples. The limit N — > oo of (2) is the Cesaro limit of the 
sequence of signs {sgn (f, K aj ,b)}j&i, which converges to af(b), if crf(b) ^ 0. This motivates the definition 
of discrete signature af on a finite set of scale samples {aj}jLi : 

af(b):=h nWb > if l^l >T < (3) 
I 0, otherwise, 

where r €]0, 1[ is an empirical threshold parameter. Interestingly, a directional statistics approach to a sign 
based signal analysis leads to the same formula (3) . This statistical derivation of the discrete signature is 
developed in the thesis of the third author [15]. 

Our discrete signature (3) has some some relations to the phase congruency of [7]. In one dimension and 
adapted to our notation, phase congruency PC is defined in [7] as the modulus of the complex quantity 

L,j=l I \J,K aj ,b) I 

Note that the normalization is carried out after summation, in contrast to our definition of discrete signature 
(2) and (3). A major issue in phase congruency is the presence of a dominating coefficient yielding a large 
phase congruency PC(6) = |PC(6)|, even though the sign values of the coefficient may point into arbitrary 
directions [7] . In [7] , this issue was resolved heuristically at the cost of introducing two additional empirical 
parameters. In our discretization (3), we do not encounter this problem since every coefficient contributes 
the same weight, no matter how small of large the modulus is. 

1.2. Organization of this article 

The structure of this article is as follows. In Section 2 we introduce the new concept of signature and 
present some of its basic properties. We then compute the signature of some specific signals and investigate 
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the behavior of signature under the fractional Laplacian and the fractional Hilbert transform. Finally, we 
set up relations to moduli based signal analyses and to the randomization of wavelet coefficients. In order 
to compute the signature for real world problems, we introduce in Section 3 a discretization of signature 
and illustrate our theoretical results by numerical experiments. 



2. Signal analysis by signature 

In this section, we introduce the new concept of complex signature for one-dimensional signals. In addi- 
tion we discuss important properties of the complex signature and present various examples. Furthermore, 
we discuss the behavior of signature under operations such as the fractional Laplacian and the fractional 
Hilbert transform and contrast it with the behavior of the classical Sobolev regularity index. We close the 
section by showing that signals, whose wavelet signs are randomly disturbed, have signature zero at every 
location. 



2.1. Definitions and basic properties 

We define the Fourier transform of a Schwartz function / G J^(1R;C) by 

&(f)(u>) :=/(w) ~Je-^f(x)dx. 

Likewise, we use the above notation for the usual extension to the space of tempered distributions S*"(W.; C). 
Furthermore, & (f) and / v denote the corresponding inverse Fourier transform of /. 
Let us introduce the class of complex wavelets we need for the definition of signature. 

Definition 1. We call a complex-valued non-zero function k G J^(R;C) a signature wavelet if k has a 
one-sided compactly supported Fourier transform, i.e., 

supp/? C [c, d], < c < d < oo, (4) 

and a non-negative frequency spectrum, i.e., 

«M > 0, for all w£l. (5) 



The wavelet system associated with a signature wavelet k is defined as the family of functions 




where a > and b £ R. 



An example of a signature wavelet is given by the inverse Fourier transform of the (one-sided) Meyer window 
W, i.e., 

k(x) = ^- 1 (W)(x). (6) 

We refer to Appendix A for the exact definition of the Meyer window function W . The graph of k is 
depicted in Figure 1. We employ this signature wavelet in the numerical experiments in Section 3. 
Recall that the sign of a complex number z £ C is defined by 



sgn z = 




if z ^ 0, 
if z = 0. 



The signature of a signal is then defined as follows. 
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Figure 1: The Meyer-type signature wavelet n (left) and its Fourier transform k (right). 



Definition 2. Let f G ,5f" (R; R) . If there exists z£C, such that for all signature wavelets n, 

limsgn(/, K a ,b) = z, 

then we define the signature, of , of f at b G R by 

a fib) := z; 

otherwise, we set 

af(b) := 0. 

Note that the signature a fib) is either equal to zero or is a complex number of modulus 1. It follows 
directly from the definition that the signature is invariant under translations, i.e., 

a{T r f)(b) = (af)(b-r) (7) 

and under dilations, i.e., 

a(D v f)(b) = iaf)iyb). (8) 
Here, the operator of translation by r G R, T r , and dilation by v G R \ {0}, D„, is defined by 

T r f(x):=f(x-r) and D u f(x):=^=f(^j, 

respectively. 

Since the Fourier transform of a signature wavelet k vanishes in a neighborhood of the origin, we have 
that 

(p, k) = 0, for any polynomial p. (9) 

Therefore, the signature is well defined on the tempered distributions modulo polynomials S/" j 2? , where 
& denotes the space of all polynomials. 

2.2. Signatures of specific signals 

In this section, we compute the signature of specific parts of a signal. We prove that the signature crfib) 
is equal to zero if / coincides with a polynomial locally around b. Furthermore, we show that the signature 
equals ±i at jumps to the right or to the left. Pure cusps yield an signature of ±1, depending on whether 
the cusp points upwards or downwards. 

Let us begin by considering some elementary examples for which the value of the signature follows directly 
from the definition. 
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Example 3. 1. Consider a polynomial p and any signature wavelet k. The vanishing of all moments of 
k implies 

(p, K a ,b) = 0, for all a > and b € K. 

Hence, polynomials have signature equal to zero on the entire real line. Note that this assertion is 
independent of the choice of k. 

2. Let us compute the signature for the cosine function: 

(cos, k 0)6 ) = (cos, (k Qi6 ) v ) = vr (5 1 + <5_i, (K a ,b) v ) = 7r v / ae" lb K(a). (10) 
Since «(a) = for sufficiently small a, we obtain 

ct(cos)O) = 0, 

for any b £ R. 

3. A similar computation applies to the sine function, and thus to all trigonometric polynomials, yielding 
the same result. 

4. In general, since the Fourier transform of a signature wavelet k vanishes in a neighborhood of the 
origin, we obtain for any band-limited tempered distribution / and for each b 6 R, 

sgn (/, K a .b) = sgn (/ v , kZJ>) = 0, for sufficiently small a. 

Hence, the signature of a band-limited tempered distribution is equal to zero for every b S K. 



To establish that the signature is equal to zero at sufficiently regular points of a signal (Theorem 5 and 
Corollary 6), we require the following lemma about the zeros of functions all of whose moments vanish. 

Lemma 4. Let g be a real-valued continuous function of rapid decay such that 

[ g(x)x k dx = 0, for all k £ N . (11) 

If the support of g is not compact, then there exists an unbounded monotone sequence {a;„} nS N of sign 
changes of g, i.e., there exists {x n }neN such that \x n \ — > oo, sgng(x2n) = 1 and sgng(x2n+i) = — L In 
particular, the set of zeros {x £ K : g(x) — 0} is unbounded. 



Proof. Assume to the contrary that there exists no such sequence. Then we can find an a > such that g 
neither changes sign on ] — oo, —a[ nor on ]a, oo[. The continuity of g implies that 



(;)"* 



) da; 



<2a max \g(x)\ =: C, Vm e N , 

x6 [—a, a] 



(12) 



where C > 0. 

Since g is continuous and not compactly supported, we find an interval [c, d] in M\ [—a, a] and a constant 
K > such that such that \g\ > K on [c, d\. Suppose then, without loss of generality, that [c, d] c]a, oo[ and 
that g > K on [c, d] . Then there exists an mo G N such that for m > toq 



.,/] 



/T\m 



dx > C. 



If 5 > on ] — oo, — a[, then choose an even m > mo, and if g < on ] — oo, — a[, then choose an odd m > mo. 
In either case, 



/ (T«<*> 



x) dx > C, 



which, together with (12), implies that 

(^-J g{x)dx>0, 

contradicting the vanishing moment condition (11). □ 

In the following, we say that a function / is of polynomial growth if there exists an Af £ N such that 
/ G Od^l^) as \x\ — > oo. Next, we show that a signal of polynomial growth has signature equal to zero at a 
point where all derivatives are equal to zero. 

Theorem 5. Let f be a real-valued, locally integrable function of polynomial growth. Further assume that 
f is smooth in a neighborhood of b S M. If 

f^{b)=0, for alike No, 

then 

af(b) = 0. 

In particular, 

supp af C supp /. 

Proof. As signature is invariant under translations (cf. (7)), it suffices to prove the claim for 6 = 0. Since / 
has polynomial growth, there exists an N 6 N such that / € 0(x N ) as \x\ — > oo. We define a function 

H(u)=u N+1 [ f(x)n(ux)dx, (13) 
Jr 

and show that H has infinitely many vanishing moments. To this end, we consider the following integrals: 
H{u)u k Au= f ( f(x)K(ux)dxu k+N+1 du = [ [ n{ux)u k+N+1 du f(x) dx, keN - (14) 



We need to justify the interchange of the order of integration. By a change of variables u — -, for x ^ 0, 
we get for the absolute value of the inner integral the following estimate: 



where Ck is a constant depending only on k. Therefore, 



\f{x)\ 



\k{ux)\ ■ M fe+iv+i du \f(x)\ dx < C k rTj^r 2 dx < oo, for all k e N . 

To justify the last inequality in (14), let e > and S > e sufficiently large. Then, since /«(0) = for all 
leN , 



Moreover, since / is locally integrable, 



6^ 



dx < oo. 



, da: < oo. 

\ x \k+N+2 
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Furthermore, as / € 0{x ) 



-6 



m d ,<oo. 



s 



Hence, the right hand side of (14) exists in absolute value and thus the interchange of the order of integration 
in (14) is justified by the Fubini-Tonelli theorem. 

Now, since all moments of the signature wavelet k vanish, we get for i^O that 



K {ux)u k+N+1 du = 0. 

JR 

Substitution into (14) yields 

/ H{u)u k du= K(ux)u k+N+1 du f(x) dx = 0, for all k G N . (15) 

JR JR J R 

We furthermore observe that H is a continuous function of rapid decay. Thus Re H and Im H are 
continuous functions of rapid decay all of whose moments vanish. First assume that H is not compactly 
supported. Then, by Lemma 4, we find an unbounded monotone sequence {w n }neN of sign changes of Re H, 
that is, 

, (+1, n even; 

sgnRci/(u„) = \ 

I —1, n odd. 

Analogously, there exists an unbounded monotone sequence {v n } n <£N of sign changes of ImiJ. Since k is a 
Hermitian function and / is real- valued, H is also Hcrmitian, i.e., 



H(-u) = H(u). 

By this symmetry property, we may assume that both sequences {u n } and {v n } are increasing and tend to 
+oo. Thus, the limit 

lim sgnH(u) 

u— >oc 

cannot exist. If H is compactly supported, then the limit equals zero. 
To conclude the proof, we observe that for u > 

sgn (Mi,„) ^( M »/«,W-,d,)^« W . 

Thus, with a = -, the limit 



lim sgn(/,/t a )= lim sgn(/, ki 

a— »U+ u— foo \ u 1 

either does not exist or equals zero. In either case, er/(0) = 0, which completes the proof. □ 

A particularly interesting consequence of Theorem 5 is the following result. 

Corollary 6. Let f be a real-valued, locally integrable function of polynomial growth which is smooth in a 
neighborhood ofb£R. If for some ko £ N , /( fe )(k) = 0, for all k>k , then 

vf(b) = 0. 

In particular, if f coincides on an open set U C R with a polynomial then o~f(b) — 0, for every b G U. 

Proof. Let p be a polynomial of order ko — 1 such that p^ k \b) — f^ k \b) for all k < ko — 1. Hence, f — p 
satisfies the hypotheses of Theorem 5 implying that a{f — p)(b) ~ 0. Since (/ — p, K a ,b) = (f,n a .b), we 
obtain er(/ — p)(b) = crf(b) = 0. The second claim is a direct consequence of the first. □ 



In the next example, we consider the unit step function. Here, we can compute the signature at b = 
explicitly. For b ^ 0, we can apply Corollary 6. 

Example 7. Let U be the unit step (or Heaviside) function defined by 



U{x) r- 

For any signature wavelet k, we have that 



1, iix>0, 
0, else. 



(U, Ka ,o) = (U, (n a ,V) = [ ^f- d^ = ^ [ de (16) 

Hence, since k > and suppK C [0, oo), we obtain that sgn (U, n a ,o) — h f° r a H a > 0. For 6^0, we apply 
Corollary 6 yielding 



i, if b = 0, 
0, else. 



aU(b) 

Note that the same arguments can be applied to the sign function x i— > 2 £/(:c) — 1. 



Next we show that the signature of a homogeneous function vanishes away from the origin. 
Theorem 8. Let f : K — > R 6e a homogeneous function of degree 7 6l. Then, 

af(b) = 0, /or any 6^0. (17) 



Proof. Let k be a signature wavelet and let 6 ^= be fixed. For any a > 0, we have that 

(f,K a ,b) = (f\Zr b ) = Va [ / v H£(a W )e-^dw = a -3 /" /V( a -i£)£(£) c -*S/«d£. 
Since / is homogeneous of degree 7, / v is homogeneous of degree —7 — 1. (See [16, Theorem 7.1.16].) Hence, 

(f,K a , b )=a^ f / v (£)£(0e-^ /a d£. 

Setting u = -, we obtain 

</, «!/„,(,) = ^ 7+ ^ / / V (0£(Oe~ <6 ** dC = «-^+^(/ v • «)(&«). (18) 

JR 

The function £ i-> k(£) ■ / v (£) is in < ^°° and vanishes on a neighborhood of the origin. Arguments analogous 
to those in the proof of Theorem 5 imply that the function u i-> ^(/ v • 1i)(bu) is of rapid decay and has 
infinitely many vanishing moments. Finally by Lemma 4, we obtain that ^"(/ v • K)(bu) has infinitely many 
zeros. Hence, 

sgn (/, K 1/u<b ) =sgn (V 7+ 5j?(/ v .«)(&u)J 

cannot converge for u — > 00 to a non-zero value. Therefore, of{b) = 0. □ 

Using Theorem 8, we may compute the signature at pure cusp-type singularities. 
Example 9. Consider the function 

f(x) = \x — xqP , where 7 > 0. 
9 



For < 7 < 1, this function resembles a pure cusp singularity. By translation invariance it is sufficient to 
compute the signature of x t— > |x| 7 . This function requires careful differentiation with respect to 7. If 7 is 
an even integer, we are in the case of polynomials and the signature is identically equal to zero. (Compare 
with (9).) If 7 is an odd integer then, according to [17, eq. 2.3(19)], we obtain 

(I • r, «.,o> = (l^R («o,o) v ) - 2(-l)^+ 1 7 ! jf {J ^W 1 d ^ (W) 
For 7 > and 7 ^ N, it follows from [17, eq. 2.3(12)] that 

(I . r, Kafi ) = (jTp, ( Ka , ) v ) = - sin (H) r( 7 + 1) de (20) 

Therefore, the wavelet signs at b = are given by 

fo, if 7 e2N , 

sgn (| . |T, KQ:0 ) = < -1, if 7 G ]0, 2[ U ]4, 6[ U . . . , 
[+1, if 7 e]2,4[U]6,8[U..., 

for every a > 0. Now, Theorem 8 implies 

cr(|.p)(6)=0, for all b ^ 0. 
Hence, by translation invariance, it follows that 

fo, if 7 e2N , 

af(x Q ) = <-l, if 7 e]0,2[U]4,6[U..., 
[+1, if 7 G]2,4[U]6,8[U.... 

and 

af{b) = 0, for b ? x . 

Having computed the signature of some basic signals, we now consider the superposition of two signals. 
Proposition 10. Let b £ M. and let /, g g ^'(R,M) be two tempered distributions such that crf(b) =t and 

(g,K ajb ) e o((f, Ko,b)), 

/or aZZ signature wavelets k. (Here, the Landau symbol o refers to a — > 0.) Then 

a(f + g)(b)=of(b). 

Proof. Let k be a signature wavelet. From crf(b) 7^ 0, it follows that (/, K a ,b) ^ for sufficiently small a. 
Therefore, 



a— >0 



K/ + s.V6>l |(/, Kaib) (i + |^) 

which completes the proof. □ 

The following example, which considers the superposition of the unit step function and a pure cusp, 
illustrates Proposition 10. 
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Example 11. Let / be the superposition of a cusp and the Heaviside function, i.e., 

f(x) = U(x) + \xp, < 7 < 1. 

From (19) and (20), we obtain 

|<|.| 7 ,Ka,0>l~ W 7+1/2 , 

and from (16), we get 

Ktf,* B ,o>|~|o| 1/2 - 

This shows that 

(I • P,K Q ,o) S o((U,K atQ )), 

for all signature wavelets k. Hence, we can apply Proposition 10 and conclude that 

af(0) = aU(0) = i. (21) 
Note that the less regular step singularity dominates over the cusp singularity. 

Next we show that, in general, a jump discontinuity induces a purely imaginary signature at the jump 
location. A function / has a jump (or step) discontinuity at b if the left-hand and the right-hand limits 
f(b+) and f(b—) exist but are not equal. 

Theorem 12. Let f be a real-valued, locally integrable function of polynomial growth and let b G R. If there 
exists a neighborhood U of b such that f is continuous on U \ {b} and has a jump discontinuity at b, then 

{-h ^f /(o-) > /(&+)• 

Proof. By the translation invariance of the signature (7), it is sufficient to prove the result for 6 = 0. First 
assume that the function / jumps upwards from the left to the right, that is, := /(&+) — fib—) > 0. 
Throughout this proof we write f s for the symmetric part and f A for the antisymmetric part of /, that 

is, 

fS{x) = /(*) + /(-«) and f A {x) = MzJk£ , 

There exists a continuous extension of f s around the origin, which we denote again by f s , so that 
lim a: _ ) .o+ / ( x ) — lima;_j.o- f S ( x ) — / S (0)- By the vanishing moments property of signature wavelets, we 
may assume that / s (0) = 0. 

Let k be a signature wavelet. Since k is real-valued, Ken a ,o is symmetric and ImK a .o is antisymmetric. 
Thus we have 

Re </,««,<)) = </ S + / A ,Re KQi0 ) = (f s ,Re Ka , ) 

and 

Im (/, n afi ) = (f s + f A ,lm K a>0 ) = </ A , Im Kafi ) . 

We observe that 

C K := / Im k(x) dx = Im / U (x)k(x) dx = Im f k(u>) dcu = [ —k(uj) duj > 0. 

Jo Jm Jr Jr u 

—= (f s ,ReK a0 ) = - / f s (x)ReK (-) dx -> 0, as a -> 0, (23) 
V a 'a J R \aJ 



We show that 
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and 



Then, 

and therefore, by Proposition 10, 



7= (f A ,lmn afi ) = - / / A (a;)ImK f-) dx ^ C K 9, as a ->■ 0. 
/a a Jj \aJ 

R-e</,Ka,o} € o(Im(/,K O)0 )) , ifa->0, 



(24) 



lim sgn (/, « O)0 ) = i. 



a->0 



To establish (23), let e > be arbitrary. Since f s (0) = and since f s is continuous around the origin, we 
can find 5 > such that \f (x)\ < e for all \x\ < 5. Thus, around the origin we get the estimate 



f s (x)ReK (-) dx <-/ Reii(-) dx = e / ° |Re«(j/)| dy < ellRe/dk!. 



Since / is of polynomial growth and Re K is of rapid decay, there exist C, M > and TV S No such that 

\f{x)\<C\x\ N and |Re«(aO| < C\x\- N ~ 2 , 
for all | a; | > M. Hence, there is an a > such that 



— M 



M 



f s (x)Re K (^) 



d;i; 



< 



C 2 



-M 



\N+2 



m) \A N+2 



dx 



< e, 



for all < a < ap. By the rapid decay property of Ren, there exists < a\ < ao, such that 



1 

- sup 

a [S,M] 



X 

Re k [ — 



- sup I Re k (x)\ < e, 

a [aS.aM] 



for all < a < a\. As / is locally integrable, it follows that 



\J-M Js J W a J-M [S,M] W 



<C'e 



with a constant C independent of a. Combining all three integral estimates, we have found that for every 
e > there exists an a\ > such that 



f(x)Re K ^) 



< (|| Rc/c|| L i + 1 + C) e, foralla<ai, 



which proves (23). 

For (24), we prove the following equivalent statement 



J ^Wlm/cQ dx-C, 



e 



0, as a — » 0. 



By the antisymmetry of Im k, this can be rewritten as 
1 



J f A (x) Im k (^ — ^ dx — J lmn(x)dx = — J (f A (x) — sgn(a;)~^ Im/t ^ — ^ da; 

Since f A (x) — sgn(z)® is continuous around the origin, (24) is proved similarly to (23). 

So far we have assumed that 9 > 0. The case 9 < 0, is treated analogously, introducing a negative sign. 
In that case, 



lim sgn (/, K a o) = -i, 

a— >0 



which completes the proof. 



□ 
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By Theorem 12, we can compute the signature of a piecewise polynomial signal, which is discontinuous 
at the interface. 
Example 13. The function 

for x > 



f{x) = 




for x < 



has a jump discontinuity to the right at zero. Thus, by Theorem 12, cr/(0) = i. Away from the origin the 
function coincides with a polynomial and, by Corollary 6, its signature is therefore equal to zero. 



2.3. Signature, fractional Laplacians and fractional Hilbert transforms 

Next, we investigate the behavior of the signature under the action of fractional powers of the Laplacian 
and the fractional Hilbert transform. We shall see that the former leaves the signature invariant whereas 
the latter acts on the signature by a rotation in the complex plane. 

We recall that fractional powers of the Laplacian (—A) 5, acting on / <E S*' (M.) / & , are defined by 

(^A)5/ := | • T • f, for rem. (25) 

We show that the signature is invariant under (—A) 5. Again, note that the signature is well defined for 
f €3?'(M.)/&>. 

Theorem 14. Let f e y'(R)/3? and r e R. Then, 

cx((-A)S/)(6)=<r/(&), forallbeR. 

Proof. Let k be a signature wavelet. By the translation invariance of the signature, we may assume that 
b = 0. Then 

((-A)i/, Ka ,„) = (|.r/,M V ) 

= a- r (/, \a . | r Ko) V ) = a- r (/, D±{\ • \ r K V )) = a~ r (/, D ((-A)5/s)> . 

Hence, as a > 0, we obtain 

sgn((-A)3/, KQ ) = sgn(/,((-A)5«) Oi0 ) . 

Since the class of signature wavelets is invariant under the application of (—A) 5, the result follows. □ 

Now we turn to the fractional Hilbert transform, which was first introduced in [18]. We follow the 
definition given in [19]. For a£K, the fractional Hilbert transform H a is defined on 5*"(R)/ £P by 

yfiTf . = e -iaf - sg n(.) . j_ (26) 

The following theorem shows that the fractional Hilbert transform T~L a acts on the signature as multiplication 
by e la i, i.e., as a rotation in the complex plane. 

Theorem 15. Let f £ y'(R)/^> and beR. Then 

a(n a f)(b) = e^ -af(b). (27) 

Proof. By the translation invariance (7), it is sufficient to prove the result for 6 = 0. Let a > 0. We have 
CH a f,^ ) = (e-^ !Sl(u »/,Ko) V ) = / e- M t s ^/( W )( Ka , ) v ( w )d W . 
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Sobolev regularity index 


Signature 


Fractional differentiation 




<r((-A)5/) = a(f) 


Fractional Hilbert transform 


s naf = Sf 


a(H a f) = e^a(f) 



Table 1: Comparison of the action of differential and convolution operators between the Sobolev regularity index and the 
signature. 



Since suppK^ c ( — °°jO]j this reduces to 

(H a f,n a , ) = [ 5 s ^)/H( KQ , ) v Hdc 

/H( M ) v Hd, = e ltt i (/,M V )=^ (/, K.,o) 



' — oo 

,15 

e 



Thus, 



a (H a f) (0) = lim sgn (H a f, K a , ) = e"** lim sgn (/, n afi ) = e M ^/(0), 

a— a— >0 



which proves the claim. □ 

A particular instant of Theorem 15 is the following result. 

Example 16. We consider the sign function f(x) = sgn(x). The Hilbert transform is given by 2 log | • |, 
regarded in the distributional sense, cf. [20, p. 80]. Theorem 15 and Example 7 yield 



a(log\.\)(b)=a(Hf)(b)=i-a(f)(b) = 



i 2 = -l, if 6 = 0, 
0, else. 



Now, let us compare the signature to the local Sobolev regularity. Recall that the local Sobolev regularity 
index s/(6) of / at b is defined by 

Sf (b) := sup {s G R : 3 ip € S>(R), ip(b) ^ 0, so that V f € W S ' 2 (R)} . (28) 

(See, for example, [13, Chapter 18.1].) The local Sobolev regularity index is often (implicitly) used in signal 
analysis since it is characterized by the decay of the moduli of wavelet coefficients, cf. [2]. Under the 
application of the fractional Laplacian, the local Sobolev regularity index of a function decreases by the 
order of the operator. More precisely, for / g 

s (-A)S/( 5 ) = s /( 6 )- r ' (29) 
see [13, Theorem 18.1.31]. On the other hand, s/ is not altered by the fractional Hilbert transform, that is, 

««-/(&) = */(&)• (30) 

To summarize, the fractional Laplacian changes only the local Sobolev index, whereas the fractional 
Hilbert transform modifies only the signature; see Table 1. In this sense, the signature can be regarded as 
being dual to classical local smoothness. 

In the examples that we have seen so far, the points of non-zero signature are a subset of the classical 
singular support. This gives rise to the question whether there is a relation between the classical singular 
support singsupp/ and the support of the signature supper/. The following two examples show that, in 
general, neither one is contained in the other. 



14 



Figure 2: The Weierstrass function for r = 0.35 and t = 9. Its singular support is the entire real line, but its signature is equal 
to for every 6 € 1. 



Example 17. Consider the WeierstraB function (see e.g. [21]) 

oo 

f(x) = r n cos(t n x), where < r < 1 and rt > 1; 



71=0 



see Figure 2. As / is nowhere differentiable, it follows that singsupp / = M. Now let k be a signature wavelet 
such that suppK c]l,£[. Then, by taking Fourier transforms, we obtain 

oo oo 

(/, K a , b ) =ttJ2 (*-t« + fa , E_ b Di_n) = v^r £ e- ibtn K(at n ). 

n=0 n=0 

Now consider the zero sequence a m — t~ m , m e N. Since supp/? c]l,t[, we have /?(i fc ) = for any k E Z. 
Therefore, 

oo 

(/, K« m ,6> = t~ % it e- w "H(t n - m ) = 0. 

Thus, if sgn (/, K a> b) converges, then it must necessarily converge to zero. Hence, 

af(b) = 0, for all bet. 

Therefore, we see that in general 

sing supp / £ supp of. (31) 



Example 18. Let / = e 1X be a Gaussian function with 7 > 0, and let k be any signature wavelet. The 
singular support of / is empty because / is smooth. On the other hand, as the support of k is not empty, 

(/,«o,o) = (/, ( K a,o) V ) = / e"^(/t Qi0 ) v M duj > 0, for all a > 0, 

implying that the signature equals 1 at b = 0. Thus, in general, 

supp of £ sing supp /. (32) 
This shows that the converse inclusion of (31) does not hold either. 



Example 18 is particularly interesting since signature considers the origin of the Gaussian as a salient 
point. This seems odd at first sight because the Gaussian function is in < ^' 00 , and thus a very regular function 
in the classical sense. However, note that the origin is the location of the global symmetry and the global 
maximum of the Gaussian, and thus may indeed be considered as a salient point. Actually, we see in the 
next section that signature can be interpreted as an indicator of symmetry or antisymmetry at a salient 
point. 
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Signal feature 


Local symmetry 


Signature 


Step to right 


antisymmetric 


+i 


Step to left 


antisymmetric 


—i 


Cusp upwards 


symmetric 


+1 


Cusp downwards 


symmetric 


-1 



Table 2: Signatures of frequently occurring feature types. We observe that the signature is imaginary at the antisymmetric 
and real at the symmetric features. 

2.4- Geometric interpretation of the signature 

Whereas the local Sobolev regularity index is by definition a measure for the local order of differentiability, 
the interpretation of signature is not so obvious. We now argue that signature can be considered as an 
indicator of local symmetry at a singular point. 

To this end, let us recall that the signature is imaginary at step discontinuities and real at cusps. This 
is a first indication that imaginary signatures correspond to locally antisymmetric features and that real 
signatures indicate locally symmetric features; see Table 2. A more rigorous argument can be derived 
from the change of the signature under the Hilbert transform. The ordinary Hilbert transform converts 
the symmetric part of a signal to an antisymmetric function, and vice-versa. Theorem 15 states that the 
application of the Hilbert transform to a signal induces a multiplication by e^" — i of its signature. Thus, 
the change of symmetry is directly reflected in the rotation of the signature by ^ in the complex plane. We 
may explicitly observe this effect for the sign function x t— > sgncc and its Hilbert transform (in a weak sense), 
the logarithm x t-¥ 2 log The antisymmetric sign function has a signature of i at the origin whereas the 
signature of the logarithm equals i ■ i = — 1. 

Finally, we notice that signature only indicates local symmetry at points with non-zero signature. How- 
ever, it is not yet clear how these points of non-zero signature can be characterized. This issue will be 
addressed in a forthcoming paper. 

2.5. Randomization of wavelet coefficients 

In this section, we show that a signal with randomly perturbed wavelet signs does not contain any salient 
points in the sense of signature. 

Let us introduce random linear operators that are diagonal in a wavelet basis. To this end, let {ipjk\j,k& 
be a compactly supported orthonormal wavelet basis of L 2 (M), where j stands for the dyadic level and k for 
the translation parameter. For / e L 2 (BLIR), we have the following decomposition 

/ = X] X ^' ^ ^ k =: X X fdkipjk, 

where convergence is in L 2 -norm. We define an operator A £ : L 2 (R) H> L 2 (R) of random perturbations of 
the wavelet signs by 

Af ■= X X tjkfjk^Pjk, (33) 
jei, fcez 

where {ejk} is a Rademacher sequence, i.e. a sequence of independently and identically distributed random 
variables taking the values ±1 with probability |. Regarding the use of Rademacher sequences in the context 
of random Fourier series, see [22]. More recently, Jaffard [23] considered the use of Rademacher sequences 
for random perturbations of wavelet series. Our next theorem shows that the signature of A £ f vanishes for 
almost every realization of the Rademacher sequence. 

Theorem 19. Let f G L 2 (K), 6 G R, and A £ an operator of random sign perturbations as defined in (33). 
Then 

a(A s f)(b) — 0, for almost every e. 
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Proof. We first compute the wavelet coefficients of A e f relative to a signature wavelet k. This yields 

sgn(A e /,K 2 -j i6 ) = sgn [^rei, m {f,ipi lW ,)(il>i m ,K2-j,b) , j > 0. (34) 
Since the sum in (33) is convergent in the £ 2 -norm for every realization of {e; m }, it follows that 

is a complex-valued bounded random variable. As {e; m } is a Rademacher sequence, we have that 

SgO-{A- e f)j,b = -sgn(AJ)^ b . 

Let I be an arc on the unit circle with positive measure. Then we have 

P(sgn(A e f) Jtb el) = P(sgn(A e f)j, b e -I), 

where P denotes the probability measure associated with the Rademacher sequence {e/ m }. This implies 
that for each zq on the unit circle 

P ( lim sgn(A e f) jtb = z b )=Q. 

Therefore, 

a(A e f)(b) = 0, with probability 1. □ 

Let us illustrate the above result by considering the characteristic function on [0,1]: / := l[o,i]- The 
signal / has jump discontinuities at and 1, but the shape of the randomly perturbed signal A e f does 
not have much in common with the original jump shape of /. However, since the amplitude of the wavelet 
coefficients is not altered, (/, ipjk) = (A e f, tpjk) > it is difficult to distinguish between / and A e f using a local 
signal analysis tool that is based only on the moduli of the wavelet coefficients. Signature, on the other 
hand, takes into account the random perturbations. The signature of the original function / is non-zero at 
and 1, but zero for the perturbed function A e f. 



3. Discretization and numerical experiments 

In practice, only a finite number of wavelet scales {a J *}^. 1 is available. Furthermore, since we cannot 
test for convergence in (1) using every signature wavelet, we have to choose a suitable signature wavelet k. 
Thus, we have to estimate the signature from the finite set of samples {sgn (/, K aj . b )}jLi- 

To motivate our numerical approach, we begin by considering the following elementary convergence result 
for discrete samples. 

Proposition 20. Let f be a tempered distribution, {oj-j-jgN a sequence such that dj — > 0, and b G R. If 
af{b) ^ 0, then 

1 N 

-^^2 s S^(f,K ai ,b) =o-f(b) (35) 
i=i 

for all signature wavelets k. 



Proof. Let k be a signature wavelet. Since crf(b) is non-zero, the sequence {sgn (/, K a ,-,b)}j&^ is convergent. 
We observe that the limit in (35) is the Cesaro sum of this sequence. The claim follows from the fact that, 
for a convergent sequence, its Cesaro sequence converges to the same limit, cf. [24, Chapter 5.7]. □ 
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Proposition 20 suggests the Cesaro limit (35) as an alternative to computing a non-zero signature af(b). 
Note that <rf(b) is of modulus 1 and so is the Cesaro limit (35). Furthermore, the elements of the Cesaro 
sequence 

1 N 

-5>gn(/,« aj , b ), iVeN, 

3=1 

are not necessarily of modulus one, but their moduli converge to 1 as N goes to infinity. This observation 
motivates the following procedure for the numerical estimation of the signature. 

Given a finite number of scale samples {fflj}^Lx' we interpret the mean of the sequence of discrete signs, 
given by 

1 N 

Wb ■= — ^sgn(/,K aj!6 ) , (36) 

3=1 

as the N-th element of a Cesaro sequence. If the absolute value |wJf,| is close to 1, we consider the Cesaro 
sequence as being convergent, with sgn Wb giving an estimate of af(b). On the other hand, a small value 
of \Wf,\ suggests a vanishing signature. More precisely, we consider \Wb\ to be non zero if it exceeds some 
empirical threshold parameter r between and 1. Hence, we propose a discrete estimate &f(b) of the 
signature of the form 

mb):=\r Wb ' * N [- T ' (37) 
10, elsewhere. 

Figure 3 and Figure 4 show numerical experiments based on the above procedure. For all experiments, 
we used the Meyer-type signature wavelet (6) and scale samples of the form aj = 2~i , with j = 0, 1, . . . , 15. 
The threshold parameter was set to t — \\f2 ~ 0.7. The first experiment (Figure 3) shows that the discrete 
estimate is in agreement with the theoretical results for pure cusps and pure steps as given in Example 9 
and Example 7. In the case of a more complex signal, we also observe that the modulus of the mean \wi,\ is 
large at the salient points; see Figure 4. These experiments show that the procedure proposed above yields 
a reasonable way to compute the signature numerically. 

Remark 1. If our sample set {sgn K aj! b)}jLi consists of non-zero elements, then the discrete signature 
(37) has some interesting connections to directional statistics [26]. In the context of directional statistics, 
the quantity Wb is called the mean resultant vector of the sample set and its sign, sgnwJ;,, the directional 
mean. Furthermore, the modulus \Wb\ of the mean resultant vector is directly connected to the directional 
variance pb via 

p h = 1 - \W b \. 

These observations allow a statistical interpretation of the discrete signature (37). The sample sequence is 
considered to be "non-convergent" if the directional variance pb is large, that is, if the orientations of the 
samples vary a lot. On the other hand, if the directional variance is small, the directional mean gives an 
estimate of the signature. 

4. Summary and outlook 

We have proposed a new tool for signal analysis that is based on complex wavelet signs to discriminate 
and ultimately classify isolated salient features. We proved that at a jump the signature is imaginary, 
whereas at a cusp it is real. Furthermore, signature is invariant under fractional Laplacians but serves as 
a multiplier when acted on by the fractional Hilbert transform. In addition, we showed that singularities 
regarded in the classical sense need not coincide with singularities described in terms of signature, and that 
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Figure 3: The discrete signature of two sample signals. We observe that the absolute value of the resultant vector is close 
to one at the salient points, here a step and a cusp, and much lower at the other points (second row). It can also be seen 
that the peaks are sharply localized, so we do not necessarily require a non-maximum suppression here. In the third row, the 
discrete signature according to (37) is depicted as a phase angle. We see that at the steps the signature is close to the angles 
i?i (right). At the cusp, the signature angles group around 7T. The experiment validates the theoretical values computed in 
Example 9 and Example 7. Note that the extra points at the boundaries are due to the periodization of the signal. 
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Figure 4: The discrete signature of a sample signal (top) taken from Wavelab [25]. We observe that the absolute value of the 
mean |«J;,| is large at the feature points and much lower at the other points (center). The bottom plot depicts the discrete 
signature in phase angle representation. We see that the signature clusters around the angles ±S at the step-like points, and 
around n and at cusp-like points. The threshold r is set equal to 0.7 in this experiment. If we choose a lower threshold 
parameter, say r = 0.4, then the discrete signature would also catch the subtle feature points, like the small jump at x = 512. 
However, in that case, we would require a non-maximum suppression to maintain the sharp localization of the pronounced 
feature points. 
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signature classifies each point of a signal whose signs are randomly distributed as non-salient. Finally, we 
proposed a method for the numerical computation and estimation of the signature for real-world signals and 
used it to validate the theoretical concepts developed in this paper. 

Future work includes the analysis of noisy signals, applications to signal reconstruction, and the multi- 
dimensional extension of the concepts introduced here. 



Appendix A. The Meyer-type signature wavelets 

A Meyer-type signature wavelet n is defined as the inverse Fourier transform of the one-sided window 
function W given by 



cos 



(ffl(5-6w)), for|<^< 



5 

3 - ~ ^ 6 ' 
5 / , , ^ 4 



W(u) := { , , ^ fOT 6<-<3, (A1) 

Vos (f g{3u - 4)) , far|<w<|, 

v else. 
Here, g is a sufficiently smooth function satisfying 

9(0=1°' Tit^ and 3(0 + 9(1-0 = 1, eeR. 
1 for £ > 1, 



For more details, see [27]. 
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